library(dplyr)
library(knitr)
library(maptools)
library(rgdal)
library(TreeSegmentation)
library(sp)
library(ggplot2)
library(rgl)
library(clue)
library(lidR)
knit_hooks$set(webgl = hook_webgl)
opts_chunk$set(warning=F,message=F)
#set color ramp for treeID
col = pastel.colors(200)
#set data paths
path_to_tiles="/Users/ben/Dropbox/Weecology/NEON/cropped_"
#set cores
cores<-3
Load in ground-truth
shps<-list.files("../data/ITCs/test/",pattern=".shp",full.names = T)
itcs<-lapply(shps,readOGR,verbose=F)
names(itcs)<-sapply(itcs,function(x){
id<-unique(x$Plot_ID)
return(id)
})
Example Pipeline
Read in Data
ground_truth<-itcs[[2]]
fname<-get_tile_filname(ground_truth)
tile<-readLAS(paste("../data/2017/Lidar/cropped_",fname,sep=""))
tile@crs<-CRS("+init=epsg:32617")
plot(tile)